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Abstract. The Drell-Yan process is a standard tool for probing the partonic 
structure of hadrons. Since the process proceeds through a quark-antiquark 
annihilation, Drell-Yan scattering possesses a unique ability to selectively probe 
sea distributions. This review examines the application of Drell-Yan scattering to 
elucidating the flavor asymmetry of the nucleon's sea and nuclear modifications to the 
CJi ' sea quark distributions in unpolarized scattering. Polarized beams and targets add 

an exciting new dimension to Drell-Yan scattering. In particular, the two initial-state 
hadrons give Drell-Yan sensitivity to chirally-odd transversity distributions. 
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O ■ 1. The Production of Massive Lepton Pairs 

Sidney Drell and Tung-Mow Yan first proposed [U Ej the process that now bears their 
names to explain a continuum of massive lepton-antilepton pairs (dileptons) that had 
been observed by Christenson et al. in proton-uranium collisions at the Brookhaven 
AGS [31 H] . The experiment was conducted to probe the large momentum transfer region 
with time-like photons to complement space-like measurements from lepton-proton deep 
inelastic scattering (DIS) data and to search for new resonances. A distinct feature of 
these data was the rapid decrease in cross section as the mass of the dilepton increased, 
as reproduced in figure [TJ While there were many plausible explanations of this spectra, 
it was the mechanism proposed by S. Drell and T.-M. Yan which described the spectra 
in terms of the (then very new) parton model of Feynman [5] that was eventually 
accepted. In this description, the dilepton cross section and its rapid decrease with 
increasing dilepton mass was explained in terms of the annihilation of a parton from 
one of the interacting hadrons with an anti-parton from the other hadron. The steeply 
falling cross section was due the paucity of large-x partons that are necessary to reach 
high mass. 
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Figure 1. The dilepton cross section as measured by Christenson et al. [J, showing a 
rapid decrease as a function of dilepton mass. The excess of events in the 3 < M 7 * < 4 
region is from the dilepton decay of the J/^. 




Figure 2. Fcynman diagram for the leading order Drell- Yan process. 



This article first reviews the basic formalism of the Drell- Yan mechanism in 
unpolarized scattering and its relation to parton distribution measurements. Next, 
angular distributions of Drell- Yan scattering and observed deviations from the expected 
distribution are reviewed. Finally longitudinally and transversely polarized Drell- Yan 
measurements are discussed. Within these discussions relevant recent and proposed 
Drell- Yan measurements will be presented. 

2. The Drell- Yan Process 

The production of massive dileptons through quark- ant iquark annihilation can be 
expressed in terms of a hard, short-distance interaction term representing the cross 
section for quark- ant iquark annihilation into virtual photon and subsequent decay to a 
dilepton pair, a qq - (illustrated in figure [2]) and the parton probability densities within 
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Figure 3. The cos 9 dependence of the proton-proton Drell- Yan cross section as 
measured by the Fcrmilab E866/NuSea experiment. The curve shows the result of a 
fit of the data to A (l + A cos 2 6) [9] . 



the interacting hadrons. The hard scattering cross section is given by 

_ 47TQ 2 1 2 

a ™~ 3M 7 2 *3 e *' [ ) 

where the cross section is reduced by the final factor of 1/3 since the color-charge of 
the quark and antiquark must match, e, is the fractional charge on the quark and M 7 * 
is the dilepton mass. To obtain the hadron-hadron cross section, it is necessary to sum 
over the available quark flavors and account for the parton distributions. To leading 
order in the strong coupling constant, a s , the Drell- Yan cross section is then 

= SrE e ' [fi(^Q 2 )fi(x2,Q 2 ) + fi(xi,Q 2 )fi(x2,Q 2 )] , (2) 

7 i 

with the sum is over quark flavors, % G {u, d, s, . . .}. The parton distributions functions 
(PDFs) are given by fi(x, Q 2 ), where x is Bjorken-x and Q 2 is the QCD scale at which 
the parton distribution is probed. In the case of Drell- Yan scattering, Q 2 = M 2 * . (In 
general, M 2 * will be used when discussing an invariant mass measured by an experiment 
and Q 2 will be used when discussing the QCD scale.) The subscripts 1 and 2 denote 
the interacting hadrons, which in a fixed target experiment, are conventionally take as 
1 for the beam hadron and 2 for the target hadron. Detailed derivations of this cross 
section may be found in the literature [SI El El- The leading order Drell- Yan mechanism 
also predicts that the spin of the virtual photon will be aligned providing a cross section 
that has a (1 + cos 2 6) dependence, where 9 is the polar angle of the lepton in the rest 
frame of the virtual photon [l], in agreement with data as shown in figure [3l Additional 
features of the angular distributions and their deviations from (1 + cos 2 6) are discussed 
in section [51 

Experimentally, one measures the momenta of the outgoing lepton and antilepton, 
allowing for the reconstruction of the virtual photon's mass, M 7 *, longitudinal 
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momentum, pi and transverse momentum, p?. It is generally more convenient to use 
the variables 

r = Ml,/s (3) 



y = i ln ( i + Pl ) ( ra P idit y)' ( 4 ) 



and 

11 = 

2 \E-p lt 

where s is the square of the center-of-mass energy of the interacting hadrons and E is 
the virtual photon's energy. From these, the momentum fractions X\ and x 2 (Bjorken-x) 
of the interacting partons are given by 

*i,2 = ( r + ^ I' 7 ' e±y (5) 



and the difference (Feynman-x) 

x F = —= « x\ - x 2 (6) 
V s 

In the limit of px — > and large \/s, this is equivalent to M 2 * = x\x 2 s and 
xf = Ipijyfs = Xi — X2- For a further discussion of the differences see [6]. 

As experiments were able to reduce the systematic uncertainty of the overall 
normalization of the Drell- Yan cross section, it became apparent that, while the leading 
order cross section (using DIS PDFs) explains many of the features of Drell- Yan 
scattering, it fails to predict the overall magnitude of the cross section by a factor 
of approximately 1.5-2. This factor, traditionally known as the "K"-factor, results from 
neglecting terms of higher order in a s in the cross section formula. The next-to-leading 
order (NLO) in a s terms [TQl IH] of the perturbative expansion, shown schematically in 
figure H] [6] , appear to account for the remainder of the experimentally measured cross 
section, within the scale uncertainty of the measurements and systematic uncertainty of 
the parton distribution input to the calculations [Hi [13], as shown in figure [5j (Recall, 
however, that while much of the global parton distribution fits is dominated by DIS 
scattering results, Drell- Yan measurement were also included in these fits.) Higher 
order QCD corrections to the cross section |T9l [20] , as well as electromagnetic radiative 
corrections to the cross section have also been calculated [2H 122] . 

The interpretation of the observed dilepton spectra in terms of parton distributions 
relies on the factorization of the Drell- Yan cross section into an infrared safe, short 
range hard scattering and the parton distributions. It further requires that these parton 
distributions have the same meaning as DIS parton distributions. In a twist expansion, 
the cross section can be expressed in terms of powers of 1/ (QR) where Q 2 is the hard 
scale and R ~ 0(1/ ^qqd) represents a non-perturbative scale [23[ 124"]. 

ff DY = ^Hard + E [l/{QR) n ] , (7) 

n 

where o"g arc j represents the convolution of the hard scattering quark- ant iquark cross 
section with the PDFs. In leading order in a s , 0"Hard ^ s §i ven by (El) but more generally 
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Figure 4. Feynman diagram for the terms of next-to-lcading order in a s for the Drell- 
Yan process: (a) and (b) QCD Compton; (d) and (e) gluon production and (f) vertex 
correction [6]. 
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Figure 5. Drell- Yan absolute cross sections measured by Fermilab E866/NuSea 
(proton-proton and proton deuterium) [I2J H3], E772 (proton-deuterium) P31 US] 
and E605 (proton-copper) [TB] compared with NLO cross section calculations based 
on the CTEQ6 [T7] and MRST [TS] parton distributions. There is an overall 6.5% 
normalization uncertainty on the E866 /NuSea data. 
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includes higher powers of a s . The sum over T n represent terms of higher twist. Collins, 
Soper and Sterman [25] 12S] and Bodwin [27] have shown that for the leading twist term, 
1/(QR)°, in this expansion the contributions of non-factorisable soft gluons cancel and 
so the leading twist term is factor isable. Qiu and Sterman showed that factorization can 
be extended to the 1/{QR) 2 term for unpolarized scattering and to 1/(QR) in a polarized 
asymmetry [23| [2lj. Factorization breaks down for terms in the expansion beyond this, 
but is not generally a problem since the Q 2 of the typical Drell- Yan experiment is greater 
than the J/ip mass squared, suppressing higher powers of 1/(QR). It may be possible, 
however, to observe the effects of power corrections by investigating regions in which 
the leading twist terms are suppressed such as high-x. Indeed, one possible explanation 
of an observed departure at large-x in pion-induced Drell- Yan from the (1 + cos 2 #) 
dependence [28l EH] [30] [31] [32] is the presence of higher twist terms as suggested by 
Berger and Brodsky [33] [33] • Drell- Yan angular distributions are discussed in greater 
detail in section [5] This explanation is consistent with the work of Qiu and Sterman [24|. 

In the t = M 2 */s — > 1 limit, the cross section is no longer adequately described 
by an expansion in a s to any fixed order. In this regime, the energies of soft and of 
collinear gluons are no longer negligible when compared with the available energy in the 
system. In the cross section, the leading-logarithmic terms of the form 

a k q - K - — -J- with (m < 2k) (8) 

s 1 - z 

are responsible for large corrections, where z = t/{x\X2) represents a partonic level 
version of r. These terms must be "resummed" to all orders in a s to adequately describe 
the process. Resummation was pioneered for the Drell- Yan process by G. Sterman [33] 
and S. Gatani and L. Trentadue [36] [37], where the resummation is described in terms of 
a Mellin transformation. A more recent alternative approach by A. Idilbi and X. Ji for 
Drell- Yan [38] [39] uses soft-collinear effective field theory, based on a similar description 
by A. Manohar [3D] for DIS, and arrives at the same result. 

The vast majority of Drell- Yan data is not near r — > 1 limit, which would 
require the interaction of very large-x parton from each hadron-a limit very difficult 
to reach when one of the partons is a sea antiquark. (Although, arguments have 
been made that resummation should be considered in any case [31], and the effects 
of neglecting resummation have been included estimates of the uncertainties in parton 
distributions [32].) Recently, the PAX collaboration [32] has proposed Drell- Yan 
measurements using an antiproton-proton collisions, thus making available valence 
antiquarks at large-x. In addition, the center of mass energy range is relatively small, 
with 30 < s < 200 GeV 2 . In these kinematics, the effects described by resummation 
will contribute substantially to the cross section [33], but are well understood. 

3. Sea Quark Distributions from Unpolarized Drell- Yan Measurements 

In a fixed target environment, where the decay leptons are boosted far forward, Drell- Yan 
scattering has a unique sensitivity to the antiquark distribution of the target hadron. 
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Combining this boost with the acceptance of the typical dipole-based spectrometer 
restricts the kinematic acceptance of the detector to xp > and consequently to 
very high values of x%, where the sea quarks are suppressed by several orders of 
magnitude compared with the valence distributions. These beam valence quarks must 
then annihilate with an antiquark in the target, thus preferentially selecting the first 
term in (J2J). This feature has been used by several recent experiments to study the 
sea quark distributions in the nucleon and in nuclei. Two new experiments have been 
proposed to extend these measurements to larger values of X2- 

Both experiments are conceptually similar to earlier fixed target Drell- Yan 
experiments. The Fermilab E906 experiment [15] will use a 120 GeV proton beam 
extracted from the Fermilab Main Injector. The experiment is already approved 
and expects to begin data collection in 2009. It will have a kinematic coverage of 
0.08 < X2 < 0.45. At the JPARC facility, a similar experiment using a 50 GeV proton 
beam has been proposed with kinematic coverage of 0.2 < x 2 ^ 0.6. An initial program 
(JPARC Phase I) using a 30 GeV beam to study J /if) physics was also proposed [46] . 



3.1. Isospin Symmetry of the Light Quark Sea 

For many years, it was believed that the proton's sea quark distributions were d-u 
symmetric, arising from approximately equal splitting of gluons into dd and uu pairs. 
While there were indications of d ^ u from Drell- Yan data in the early 1980 's [17] 
it was observation of a violation of the Gottfried Sum Rule [IE] in muon DIS by the 
New Muon Collaborations [49j [50] that forced this belief to be reconsidered. Because 
of its sensitivity to the antiquark distributions Ellis and Stirling suggested [51] using 
Drell- Yan as a probe of the light quark flavor asymmetry with hydrogen and deuterium 
targets. To illustrate this sensitivity, in leading order with the X\ 3> X2, assuming charge 
symmetry and the dominance of the uu annihilation term, the ratio of the per nucleon 
proton-proton to proton-deuterium Drell- Yan yields can be expressed as 

The NLO terms in the cross section provide a small correction to this ratio and are 
considered in the analysis of the data, along with the deviation from the X\ ^> X2 
limit. Fermilab E772 used their existing proton induced Drell- Yan data to compare 
the W/C and W/ 2 H yields to extract upper limits on the isospin asymmetry for 
0.04 < x < 0.27 [52]. The first dedicated measurement of the d/u asymmetry using 
Drell- Yan scattering was made by the CERN NA51 experiment [53]. The acceptance 
of the NA51 toroid-based detector was such that the average rapidity (y) = and 
x\ = X2 = 0.18. The asymmetry measured by NA51 was 

pp a pp _ pn 

a by = 2 7^ - 1 = aPP + aPn = ~ om ± °- 02 ( stat -) ± °- 025 ( s y st -)> ( 10 ) 



0~pd 



2(7 pp 



2 
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= 0.51 ± 0.04 (stat.) ± 0.05 (syst.), (11) 

i=0.18 



with the second equality only valid if nuclear effects are ignored. From this NA51 
extracted 

u 
d 

a clear signal for isospin symmetry violation in the sea antiquark distributions. 

The Fermilab E866/NuSea experiment used Drell- Yan to measure x-dependence of 
the dju ratio. Fermilab E866/NuSea used a spectrometer composed of three dipole 
magnets. The first two magnets served to focus large transverse momentum, p?, 
dimuons into the spectrometer while tracking surrounding the third magnet provided 
a momentum measurement of the individual muons. The experiment used 800 GeV 
protons extracted from the Fermilab Tevatron incident on hydrogen and deuterium 
targets. The remainder of the beam which did not interact in the targets was intercepted 
by a copper beam dump contained within the first magnet. Additionally, the entire 
aperture of the first dipole was filled with copper, carbon and borated polyethylene, 
absorbing essentially all particles other than muons produced in the interaction of the 
beam with the targets or beam dump. (Fermilab E772 p3] used essentially the same 
apparatus.) E866/NuSea recorded 360,000 Drell- Yan events, approximately two thirds 
from a deuterium target and the remainder from a hydrogen target. The ratio of Drell- 
Yan cross sections, a pd j (2a pp ), measured by E866/NuSea as well as the extracted ratio 
d(x)ju(x) is shown in figure E] When these cross section ratios were included in global 
parton distribution fits [HI [551 EE] , they completely changed the perception of the sea 
quark distributions in the nucleon. 

The E866/NuSea data present an interesting picture of the sea quark distributions 
of the nucleon that may shed some light on the origins of the sea quarks. At moderate 
values of x the data show greater than 60% excess of d over u, but as x grows larger, this 
excess disappears and the sea appears to be symmetric again. If the sea's origins are 
purely perturbative, then it is expected to have only a very small asymmetry between d 
and u [60| IB"!] . Many non-perturbative explanations for the origin of the sea including 
meson cloud models [621 [63] , chiral perturbation theory [611 ES] or instantons [661 EI] 
have been suggested which can explain a large asymmetry, but not the return to a 
symmetric sea which is seen as x — > 0.3. These models are reviewed elsewhere [68] • It is 
interesting to note the importance of dynamical chiral symmetry breaking that Thomas 
et al. have related to the d — u difference through the presence of Goldstone bosons 
which form a pion cloud around the nucleon [69]. Henley et al. have calculated the 
coordinate-space distribution of the d — u asymmetry and observe that this distribution 
agrees with what is expected if d — u originates from a pion cloud surrounding the 
nucleon [70]. Plotted in terms of d(x) —u(x), the observed E866/NuSea asymmetry, as 
shown in figure [7] [57] , can be compared directly to non-perturbative models since the 
(flavor symmetric) perturbative component of the sea is removed. 

The E866/NuSea data become less precise as x increases beyond 0.25 and the 
exact trend of dju is not clear. To help understand this region better, Fermilab E906 
experiment has been approved to collect Drell- Yan data in this region. The Fermilab 
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Figure 6. (a) The blue squares show the ratio the proton-deuterium to twice 
the proton-proton Drell- Yan cross sections versus x-i as measured by Fermilab 
E866/NuSea [57] . (b) The blue squares show d(x)/u(x) ratio extracted by 
E866/NuSea [57] The magenta triangle is the NA51 [53] measurement of d/u. The 
central curve in the cross filled band shows the d/u ratio from the CTEQ6m fit [17], 
which included the E866/NuSca data, and the band represents the uncertainty from 
the fit. The curves labeled CTEQ4M [58] and MRS(r2) [59] show the parameterizations 
of d{x)/u(x) which included the NA51 point and the EMC integral but not the 
E866/NuSea data. 




Figure 7. The d(x) — u{x) distribution as extracted by E866/NuSea (blue squares) 
using Drell- Yan [57] and by HERMES (green triangles) using semi-inclusive DIS [71] . 
Also shown are a pion model calculation of Peng et al. [53] based on the procedure of 
Kumano [55] , 
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E906/Drell-Yan experiment [15] is modeled after its predecessors, Fermilab E772 and 
E866/NuSea. The E906 experiment will use a 120 GeV proton beam rather than the 
800 GeV beam used by E866 and E772. Experimentally, the lower beam energy has 
two significant advantages. First, the primary background in the experiment comes 
from J/ip decays, the cross section of which scales roughly with s. The lower beam 
energy implies less background rate in the spectrometer and allows for a correspondingly 
higher instantaneous luminosity. Second, the Drell-Yan cross section at fixed X\ and x^ 
is inversely proportional to s (recall pa x\X2s) and thus the lower beam energy 
provides a larger cross section. The muons produced in a 120 GeV collision have a 
significantly smaller boost, which forces the apparatus to be shortened considerably in 
order to maintain the same acceptance. The smaller boost will also create a larger 
background from pion decay. Fermilab E906 expects to have statistical precision better 
than 1% for x < 0.35 and 10% for 0.35 < x < 0.45, a clear improvement over the 
E866/NuSea data. 

A similar experiment jl6] has also been proposed for the JPARC facility. This 
experiment would employ a 50 GeV proton beam with an apparatus similar in design 
to the E906/Drell-Yan apparatus. The initial phase of the JPARC facility is for a 30 
GeV synchrotron, which kinematicly has insufficient phase space above the J/ip in mass 
for a Drell-Yan experiment; although, once the entire facility is completed, including 
capabilities for 50 GeV beam, a significant Drell-Yan program could be mounted. 

3.2. Antiquark Distributions of Nuclei 

The distributions of partons within a free nucleon differ from those of a nucleon bound 
within a heavy nucleus, an effect first discovered by the European Muon Collaboration 
(EMC) in 1983 [72]. These nuclear effects are now generally divided into four regions 
in x-space: the shadowing region with x < 0.1, the anti-shadowing region covering 
0-1 ^ x ^ 0-3, the EMC effect region for 0.3 < x < 0.6 and finally a region dominated 
by Fermi motion of the nucleons with 0.6 < x. (For a review of the nuclear EMC 
effect, see [73].) Almost all of the data on nuclear dependencies is from charged lepton 
DIS experiments, that are sensitive only to the charge-weighted sum of all quark and 
antiquark distributions. Nuclear effects in the sea quark distributions may be entirely 
different from those in the valence sector |74| . but an electron or muon DIS experiment 
would not be sensitive to this. With the ability to probe the antiquark distributions of 
the target, Drell-Yan presents an ideal tool with which to distinguish between sea and 
valence effects. 

Some early models of the EMC effect were based on the convolution of a virtual pion 
cloud with the nucleon. In these models, virtual pion contributions to nuclear structure 
functions were expected to lead to sizable increases in sea distributions of the nuclei 
compared with deuterium [751 C23 EZ]- Early Drell-Yan studies at both CERN [78j EH] 
and Fermilab [17] lacked the statistical sensitivity to observe the expected effects. In 
1990, measurements by Fermilab E772 [51] found that the expected enhancement was 
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Figure 8. E772 measurements of the ratio of Drell- Yan cross sections on (a) carbon, 
(b) calcium, (c) iron and (d) tungsten to deuterium [54] . Aside from shadowing in 
the X2 < 0.2 region, no nuclear effects are observed. In (c), to illustrate the level 
of effects expected, curves based on several different representative models are also 

plotted [751 m [ID [5H [H[E3 [El- The differcnccs betwcen these curves is discussed 
briefly in the text. 



clearly absent, as shown in figure [HI The non-observation of evidence of a pion excess 
calls into question the most widely believed traditional meson-exchange model [80] of 
the nucleus. The expected enhancement to the sea is illustrated in Fig. [8b, which shows 
the expected Drell- Yan ratio in iron to deuterium, based on nuclear convolution model 
calculations by Coester [75J [76], [77] meant to explain the originally observed EMC effect. 

More recent calculations [EH E2] EH EI], made in light of the E772 data, predict a 
smaller nuclear dependence, consistent with the statistical uncertainties of E772. Jung 
and Miller [HI] revisited the calculations of Berger and Coester and examined the effect of 
the quantization of the pions on the light cone versus at "equal time" . With "equal time" 
quantization, they calculate a roughly flat 8% increase in the Drell- Yan iron cross section 
over the deuterium cross section. Brown et al. [H2] argue that with the partial restoration 
of chiral symmetry, the masses of hadrons made up of light quarks decrease with density. 
This rescaling leads to altered couplings, which lead to an overall x-dependent decrease 
in the Drell- Yan cross section in nuclei. Dieperink and Korpa [HH] also argue that the 
Drell- Yan cross section ratio should decrease in a nucleus. Their arguments are based 
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on particle- and delta-hole model, which results in a strong distortion of the free pion 
structure function. Based on the chiral quark soliton model, Smith and Miller predict 
essentially no nuclear dependence for Drell-Yan, while at the same time qualitatively 
explaining the observed EMC effect [S3]. For x > 0.2, the E772 statistical uncertainties 
allow some freedom for these models. At x ~ 0.3 these newer models have nuclear 
effects of the order 5 to 15% in the Drell-Yan ratio and tend to diverge from each 
other. Fermilab E906 will provide the sensitivity needed to see these differences with 
an expected statistical precision of 1% at x = 0.3 and data out to x — 0.45 [32]. 

4. Unpolarized Parton Distributions at Large- x of the Beam Hadron 

In addition to probing the sea quark distribution of the target nuclei, Drell-Yan 
scattering can be used to probe the structure of the beam hadron. As shown in (T5]), the 
cross section is a convolution of the beam and target parton distributions. The same 
feature of the acceptance which allows fixed target experiments to probe the target's sea 
quark distributions gives experimental access to the high- a;, valence parton distributions 
of the incoming hadron. Absolute cross section measurements have been used to explore 
both pionic and protonic parton distributions; although, some authors suggest that 
there may be breakdown in factorization at large-x [33l 134"] . Such effects could cloud 
the partonic interpretation of these data. 

The large-x parton distributions of the proton are relatively unknown, both in 
absolute magnitude and in the ratio of d/u. Experimentally, these have been accessed 
through DIS. Unfortunately, much of the high-x data used in DIS measurements 
involve nuclear targets with relatively unknown and possibly large nuclear corrections 
(even to deuterium) [85]. Proton induced Drell-Yan, which is sensitive to the 4u + d 
distribution of the beam proton, is an alternative way to reach high-x with no nuclear 
corrections. The proton-proton and proton-deuterium cross section measurements 
of Fermilab E866/NuSea when compared with NLO calculations based on both the 
CTEQ6m [17] and MRST [55] parton distributions show a small but systematic trend 
suggesting an overestimation of the strength of this parton distribution at large-x 
(x < 0.8 however) [12]. Future Drell-Yan experiments such as Fermilab E906 [35] 
or a JPARC-based experiment [36] will be able to improve dramatically on these 
measurements in both statistical precision and reach in x. 

The pion's parton distributions are of considerable interest because of the pion's 
many unique roles in nuclear physics. The pion's valence antiquarks have been used to 
explain partially the observed d/u ratio in the proton through a pionic cloud around 
the bare nucleon. (See Sec. 13.11 ) Models of nuclear binding depend on pion exchange, 
and so the pionic quark distributions should modify the parton distributions in nuclei. 
(See Sec. 13.21 ) While the pion is a qq system, its extremely low mass arises from its role 
as the Goldstone boson of dynamical chiral symmetry breaking and this role must be 
considered in any descriptions of the pion's partonic nature. 

There are a variety of theoretical and model-based descriptions of the pion's valence 
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parton structure. At large-x valence structure of the pion can be parameterized as 
(1 — Xf, where QCD evolution of the parton distributions makes (3 a function of Q 2 . At 
low Q 2 , arguments based on the parton model [86], pQCD [87J, EE] and Dyson- Schwinger 
equations [89] require that /3 ~ 2. At the same time the Drell- Yan- West relation [90l [9Tj . 
duality arguments [92] and Nambu-Jona-Lasinio models [93j EU [95], [96] favor a linear 
dependence of ~ 1. An early leading order analysis of pion induced Drell- Yan data 
found (3 = 1.26 ± 0.04 [32] with (M r ) = 5.2 GeV. Noting that the strength of higher 
order diagrams could have a kinematic dependence, a fit to the same data, this time in 
NLO found (3 = 1.54 ± 0.08 [97], still somewhat in between a linear and quadratic 1 — x 
behavior. Some future pionic Drell- Yan experiment with improved kinematic resolution 
in the high-xi region could help to resolve this question. 

5. Unpolarized Angular Distributions 

In leading order, ignoring transverse momenta (kx) of the interacting partons, the 
Drell- Yan angular distribution is naively expected to have the form (1 + cos 2 #). More 
generally, Collins and Soper [98] have shown that the expression for the angular 
distribution is 

— oc 1 + A cos 2 9 + /i sin 29 cos H — sin 2 9 cos 20, (12) 
dvi 2 

where 9 is the polar angle of the positive lepton in the rest frame of the virtual photon 
and is the azimuthal angle. The additional terms arise from the &t of the interacting 
partons and higher order graphs in a s . After consideration of the intrinsic kx of the 
partons, care must be taken in precisely defining the the z-axis of rest frame of the virtual 
photon. The most common choices for this definition are the w-channel frame where 
the the z-axis points anti-parallel to the target nucleon direction; the Gottfried- Jackson 
frame (t-channel) [99j has the z-axis pointing parallel to the beam nucleon and the 
Collins-Soper [98] frame where the z-axis bisects the angle between the w-channel and 
Gottfried- Jackson z-axes, in an attempt to minimize the effects of kx on the observed 
angular distributions. The transformation between A, fi and v in the three frames is 
straight-forward [3Tj . 

In NLO, a relationship between A and v of the general angular distribution formula 
in ( fT2"l) was derived by C.S. Lam and W.-K. Tung [lOOj . In analogy to the Callan-Gross 
relationship of DIS |101j , the Lam- Tung relation states that 

1 - A = 2za (13) 

Unlike the Callan-Gross relation, the Lam- Tung relation is expected to be largely 
unaffected by QCD [100] . including resummation effects [102j . 

The validity of the Lam- Tung relation has been studied with both pion-induced 
Drell- Yan by CERN NA10 [301 [29] and Fermilab E615 [3H [32] and proton induced 
Drell- Yan by Fermilab E866/NuSea [103] . Pionic Drell- Yan experiments have observed 
a violation of the Lam- Tung relation. This violation is most prominent at high transverse 
momentum of the dilepton, px, where v rises without a corresponding decrease in A, as 
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Figure 9. The parameters (a) A, (b) /j, and (c) v [see (fT3|) ] of fits to the angular 
distributions as a function of pt from pion (CERN NA10 [30} US] and Fermilab 
E615 [HDIH2]) and proton (Fermilab E866 [103]) induced Drell- Yan scattering. For 
the NA10 experiment, the data set with the best statistical precision, 194 GeV, is 
shown. Plotted in (d) is 2v — (1 — A) which is predicted to be zero by the Lam- Tung 
relation. For the pionic Drell- Yan data there is evidence for a violation of the Lam- 
Tung sum rule, especially at large-p^. The uncertainties are statistical only and in (d) 
do not include possible correlations between A and v. 



shown in figure [9] The violation appears to be independent of the target nucleus [29J . 
In contrast, the recently released Fermilab E866 proton induced Drell- Yan data are 
consistent with the Lam- Tung relation, even at high-pr [103J. 

Berger and Brodsky have suggested that as x — > 1, higher-twist effects cause 
the polarization of the virtual photon to change from (1 -l-cos 2 ^) to a ^-dependent 
sin 2 9 [33] . While effects consistent with this have been observed in pionic Drell- Yan 
data [32] at high-x, the violation of the Lam- Tung relation is seen over a much broader 
range in x than would be expected by this explanation. Brandenburg, Nachtmann and 
Mirke hypothesized a spin correlation in the violation of factorization that would give 
rise to a nonzero cos 20 distribution [104] . which could explain the pionic Drell- Yan 
data. 

An alternative explanation proposed by Boer [105j based on the existence of a chiral- 
odd, T-odd distribution function, h^(x,kT), with an intrinsic transverse momentum 
dependence, k?, of Boer and Mulders [106] . This distribution function is an analog of 
the Collins fragmentation function [107] . It represents the correlation of the parton's 
transverse spin and k^ in an unpolarized nucelon. Boer argues that the presence of 
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the hi(x,kT) distribution function will induce a cos 20 dependence to the Drell-Yan 
cross section and fits the observed NA10 [30] data to a crude model of this distribution 
function. Within a quark spectator-antiquark model, Lu and Ma have shown that 
the observed cos 20 distribution can be reproduced with nonzero hj;(x, kr) in both the 
pion and target nucleon |108j . Other transversity distributions are discussed again in 
section [61 

When considering any of these explanations for the violation of the Lam- Tung 
relation in pion-induced Drell-Yan, it is important to remember that these results must 
be reconciled with the apparent absence of a violation in proton-induced Drell-Yan. One 
significant difference is that the valence antiquark in the pion-induced case allows the 
experiment to probe the quark distributions of the target, while in the proton induced 
case, only the target antiquark distributions are probed. Alternatively, the possible 
interpretation as a higher-twist effect might have a kinematic dependence on y/s [104] . 
Recall that the pionic data had s/s = 11 and 16 GeV while the protonic data had 
\fs = 39 GeV. Such an effect should clearly be seen then in the upcoming Fermilab 
E906 experiment with proton-induced Drell-Yan at -^/i = 15 GeV. 

6. Transversity Measurements with Drell-Yan Scattering 

Polarized beams and targets add an extra dimension to Drell-Yan scattering, enabling 
it to be used to study both quark helicity (longitudinal) and transversity distributions. 
Close and Sivers suggested that the sea quarks produced through gluon splitting 
would be polarized and that this effect could be observed with longitudinally polarized 
Drell-Yan scattering |109j . Ralston and Soper considered the possibility not only of 
longitudinally polarized but also of transversely polarized scattering, proposing the 
measurement of certain asymmetries as a test of the Drell-Yan Model [110] . 

In the near future, there are several proposed experiments at polarized beam 
facilities. One of the options that may be available in Phase II of the JPARC program 
would include a polarized proton beam. The FAIR program [lllj at GSI includes 
plans for a polarized antiproton ring. The PAX collaboration has proposed [13] a 
series of Drell-Yan and charm measurements using this beam. The first phase of the 
experiment would involve polarized (or unpolarized) antiproton beams of up to 3.5 GeV 
colliding with an internal polarized hydrogen target. In the second phase, the PAX 
experiment will be run as a collider experiment with up to 15 GeV polarized antiprotons 
colliding with a beam of 3.5 GeV polarized protons. The RHIC-spin program also 
offers opportunities for measuring polarized Drell-Yan scattering either with the existing 
detectors or a new detector dedicated to polarized Drell-Yan physics [1121 1113j . 

In leading twist, in addition to the unpolarized parton distributions, /, and the 
helicity distributions, Af, a complete description of the proton requires knowledge 
of the transversely polarized parton distributions, denoted h\. (Various notations for 
transversity distributions are used in the literature. This notation corresponds to that 
of Jaffe, Ji and Mulders [1141 1115] . Alternatively, A^/ in the notation of Barone 
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and Ratcliffe [1161 1117] or 5q have been used.) The quantity hi represents the net 
transverse helicity of the quarks in a transversely polarized hadron. This distribution is 
just as fundamental as the unpolarized and helicity distributions; although, much less 
studied. While the existence of the transverse parton distributions was recognized in 
the late 1970's [US] , they received little consideration in discussions of proton structure 
until recently. This is, perhaps, largely because of the difficulty in measuring these 
distributions. Because transversity is a chirally odd quantity, it cannot be probed 
in inclusive DIS, the traditional tool for deducing parton distributions. In order to 
observe hi chirality must be flipped twice |118j . This requires two hadrons in the 
interaction, for example, one each in the initial and final state or two in the initial state. 
Thus, both semi-inclusive DIS and Drell- Yan scattering offer excellent opportunities to 
study transversity distributions. In the spin-| nucleon, there is no gluon transversity 
distribution and so the QCD evolution of hi is quite different from the QCD evolution 
of Af. The distribution hi evolves as a flavor non-singlet quantity [117] . 

With the additional consideration of intrinsic transverse momentum, kx, in the 
proton, two other distributions emerge. The Boer-Mulders function, /^(cc, /cr) [106] . 
was introduced in section [5] to explain the observed cos 20 distributions in Drell- Yan 
scattering. The Sivers distributions function, f^p(x, £>r), characterized a correlation in 
a polarized nucleon of the unpolarized parton density with kj< |119l 1120] , and can be 
observed in single spin asymmetries as discussed in section I6TT1 While these distributions 
functions vanish when integrated over kx, their possible existence and /^-dependence 
offers a window into spin-orbit correlations in the nucleon. Recent measurements of 
semi-inclusive DIS have provided an initial insight into both transversity distributions 
and fragmentation functions; although disentangling the effects of fragmentation from 
the distribution functions can be difficult. With Drell- Yan scattering, the fragmentation 
functions are not relevant and one has a clean probe of only the distribution functions. 
For excellent reviews of the physics of transversity, see [1161 1117] . 

6.1. Single Transverse Spin Asymmetries and the Sivers Function 

In hadron-hadron interactions with one of the two hadrons is transversely polarized and 
the other is unpolarized, surprisingly large asymmetries of the form 

^ _ da^ - da 1 
da^ + da* 1 ' 

have been observed {e.g the Fermilab E704 experiment [121] which measured p^p — > 
irX). D. Sivers suggested (prior to the publications of the E704 results) that these 
single transverse-spin asymmetries (SSA) could be used to provide information on the 
dependence of unpolarized partons in a transversely polarized nucleon and that 
these SSA were non-vanishing when the transverse motion of the individual partons 
is considered [1191 ri20j. He applied these argument to the reaction hp^ —> irX, where h 
represents any unpolarized hadron, p^ is a transversely polarized target and ir represents 
any detected spinless final state meson. The asymmetry arising from unpolarized 
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quarks in a transversely polarized hadron is known as the "Sivers asymmetry" with the 
associated "naively" T-odd f^(x,kT) parton distribution. Initially, it was incorrectly 
believed that An would vanish because of time- reversal invariance [107] . In QCD 
at leading twist it is also expected to be zero [122] . It was later shown that with 
correct consideration of Wilson lines that the "Sivers asymmetry" can be non- vanishing 
at leading twist. The Wilson lines provide for the gauge invariance in the parton 
number densities and appear as soft initial-state interactions in Drell- Yan [1231 1124] . 
Alternatively, with collinear QCD factorization at twist-three, it was recognized that 
there could be a non-zero SSA [1251 1126"] . Both mechanisms, the fey-dependent parton 
distribution approach and twist-three QCD approach, provide mechanisms to explain 
the same physical observable-a large SSA. Within the context of Drell- Yan scattering, 
the fey-dependent parton distributions are best applied to the domain in which px *C M 7 * 
while the twist-three QCD approach is more applicable for A QCD <C px- Recent work 
by Ji et al. has shown that these two approaches are in fact related, and that they give 
the same result in the overlap region, A QCD <C pr *C M 7 * [1271 1128] . 

The single spin asymmetries observed in the Fermilab E704 experiment [121] could 
either be attributed to the f£r(x, fey) distribution in the proton or to a transversity 
dependent "Collins" fragmentation function [107] . While these two mechanism have 
similar signatures in the E704 experiment, they produce different angular distributions 
in semi-inclusive DIS. Recently, the first experimental evidence for both the Sivers 
f^p{x,krp) distribution function and the Collins fragmentation function has been 
observed by the HERMES experiment [129] . Both the Sivers distribution and Collins 
fragmentation functions have been extracted with global analyses [1301 1131] based on 
observations from HERMES p9] and COMPASS [132]. These analyses and model 
calculations [133] predict asymmetries of approximately 8% in p^p Drell- Yan in the 
kinematics of the proposed PAX experiment [13] , well within their statistical precision. 
For p^p Drell- Yan at RHIC, these analyses predict a 1-10% asymmetry on the kinematics 
accepted. One very interesting consequence is that the Sivers asymmetry should have 
the opposite sign when measured in Drell- Yan scattering [123] . i.e. 

f$r(x, k T )\ ms = - fi T (x, k T )\ DY ■ (15) 

The experimental verification of ([TBI) is a key to validating the current understanding 
of kx effects in transversity distributions and is one of the goals of the next generation 
of transversely polarized Drell- Yan experiments. 



6.2. Double Spin Asymmetries 

Doubly, transversely polarized Drell- Yan scattering offers the possibility of measuring 
hi(x) without any complications from fragmentation. The asymmetry is given 

by [nm ma msi 

da^ - da n 
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J2i e l [hii(xi)h u (x2) +hu(xi)hu(x2)] , . 

where 

sin 2 6 cos (20) . . 

«tt = — — — fir- 17 

1 + cos" 1 # 

Note that this represents a convolution of the quark hi(x\) and antiquark hi{x 2 ) 
transversity distributions. Unlike unpolarized Drell-Yan scattering, which also measured 
a quark-antiquark convolution, the assumption that the quark distributions are known 
from previous data (primarily DIS) is no longer valid and so a significant amount of 
data is required to untangle this convolution. Furthermore, because the gluons do not 
directly couple to h\, it is expected that the sea contributions to the asymmetry may 
be small. 

As an alternative to this, the PAX collaboration [42] has proposed the measurement 
of antiproton-proton Drell-Yan scattering in a collider mode. In this case the cross 
section is dominated by valence antiquarks annihilating with valence quarks. Assuming 
the dominance of the w-quark terms, the asymmetry in (Tl6l) reduces to |43j 

a dy „ hi u (xi)hi u {x 2 ) , . 

A TT = a TT — — — — , (18) 

u{xi)u{x 2 ) 

thus giving a direct experimental access to \hi u \. Statistically, this is a challenging 
measurement, as polarized antiprotons are difficult to create. By relaxing the usual 
mass requirement for Drell-Yan to include events both above and below the J/ip, PAX 
hopes to be able to obtain statistical precision of 10%. The relatively low energy of PAX 
places some data in a region in which resummation may be important in interpreting 
the overall cross sections. Fortunately, it appears that the asymmetry is largely 
unaffected by resummation 



7. Measurement of Quark Helicity with Longitudinally Polarized Drell-Yan 

The total spin of the proton can be decomposed into the spin and orbital contributions 
from the quarks and gluons: 

(s?) = l = l*Z + L q + AG + L gi (19) 

where AS and AG represent the contributions of the quark and gluon helicity 
respectively and L q and L g are the orbital angular momentum of the quarks and gluons. 
Motivated by the observation of the European Muon Collaboration (EMC) that only a 
very small fraction of the total spin of the proton is carried by the quarks [1341 1135j . 
an analysis of the world's data by the SMC group finds the fraction of the proton's 
spin carried by quarks at AS = 0.38 ± 0.06 [136J. The quark spin, AS, can further be 
decomposed into contributions of the individual flavors of quarks and antiquarks. 

The HERMES experiment at DESY has been extremely successful at studying 
the flavor decomposition of AS using semi- inclusive DIS |137j . The results of this 
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type of experiment depend on precise knowledge of both the polarized and unpolarized 
fragmentation functions. Some authors have argued that the contributions of the sea 
quarks to the polarized asymmetries are small, reducing the sensitivity of this type of 
experiment |138[ 1139] . In Drell- Yan scattering, however, this is not the case [138[ 1140] . 
and knowledge of the fragmentation functions is not necessary. 

As with unpolarized Drell- Yan, polarized Drell- Yan scattering offers a window into 
the sea quark distributions. The asymmetry for longitudinally polarized Drell- Yan, A^ L 
is given by [109 J 

da ++ - da + ~ 



da++ - da+- 

E,e? [Af i (x 1 )Af i (x 2 ) + A/i^OA/i^)] 



(20) 



Ei e f [fi(xi)fi(x 2 ) + /i(xi)/i(z 2 )] 
Here, da ++ (da + ~) denotes the spin parallel (anti-parallel) Drell- Yan cross sections. 
Making the same "fixed target" assumptions as in the unpolarized case, ( J20l) shows that 
measures Am in the target proton, convoluted with the large- £ Au distribution in 
the proton. In the limit of exact SU(6), Close and Sivers [109] have shown that 

^K^H^^f (21) 

With the addition of a polarized deuterium target, it would be possible to extract Ad 
as well, providing that the valence spin distributions are known. Such a measurement 
has already been proposed for the J-PARC Phase II facility with a 50 GeV polarized 
proton beam. The proposed experiment expects to be able to achieve 10% statistical 
precision for 0.3 < x < 0.5 [46 J. 



8. Conclusions 



Over the last 35 years, Drell- Yan scattering has played an important role in elucidating 
the hadronic structure and will continue to do so with the arrival of several new polarized 
and unpolarized experiments. Unpolarized Drell- Yan scattering has been critical in 
measuring the sea quark structure of the nucleon and the nucleus. Fermilab E866 /NuSea 
has probed the flavor asymmetry of the antiquark sea, showing conclusively that non- 
perturbative processes contribute to the sea. Fermilab E772 measured the nuclear 
effects on the sea quark distributions, finding that widely accepted models of nuclear 
binding to be lacking. With transversely polarized beams, these experiments will enable 
the extraction of the hi(x) transversity structure function. In combination with DIS 
scattering, transverse single spin asymmetries will test the understanding of the kr 
dependent Sivers function, f£p(x, Ut)- 
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